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SUMMARY 


This report describes the latest addition to the STARS system of computer programs, STARS-2P, 
for the plastic, large deflection analysis of axisymmetrically loaded shells of revolution. The STARS 
system uses a numerical integration scheme to solve the governing differential equations. Several 
unique features for shell of revolution programs that are included in the STARS-2P program are de- 
scribed. These include orthotropic nonlinear kinematic hardening theory, a variety of shell wall cross 
sections and discrete ring stiffeners, cyclic and nonproportional mechanical and thermal loading cap- . 
ability, the coupled axisymmetric large deflection elasto-plastic torsion problem, an extensive restart 
option, arbitrary branching capability, and the provision for the inelastic treatment of smeared stiff- 
eners, isogrid, and waffle wall constructions. To affirm the validity of the results, comparisons with 
available theoretical and experimental data are presented. 
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LIST OF SYMBOLS 


Lower Case Latin: 

e n Linear strain 

{f} Local force matrix 

h Shell thickness 

k Shell curvature 

[k] Stiffness matrix 

m Distributed moment load 

n Fourier harmonic number 

0 Subscript meaning reference surface 

r t Meridional radius of curvature 

r 2 Circumferential radius of curvature 

r c Ring centroidal radius 

r 0 Radius of revolution 

s Arc length 

u Circumferential shell displacement 

v Meridional shell displacement 

w Normal shell displacement (positive inward) 

Upper Case Latin: 

C Eccentricity of reinforcement 

D Bending stiffness 

E Young’s modulus 

F Distributed loads on the shell 

G Shear modulus 

J Effective transverse shear stress resultant 

K Extensional stiffness 

M Bending stress resultant 

N Membrane stress resultant 

Q Transverse shear stress resultant 

T Temperature 

Z, R, 6 Global coordinate system 
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Greek 


a Coefficient of thermal expansion 

w Rotation of the normal to the middle surface of the shell 

£ Shell coordinate normal to the middle surface of the shell (positive inward) 

6,<p Circumferential and meridional shell coordinates 

Oj 5 Normal components of stress 

v Poisson’s ratio 

Tj j Shearing components of stress 

A Deformation 

{6} Local displacement matrix 

Other symbols are defined in the text where used. 
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INTRODUCTION TO NUMERICAL INELASTIC SHELL ANALYSIS 


Recent innovations in digital computer technology allow designers to analyze shell structures of 
complex configurations without unduly restrictive approximations. A number of versatile computer 
programs, based on various methods of analysis, are presently available for the inelastic analysis of 
shells of revolution. These are briefly discussed, and categorized by the numerical method used, 
below. 

Finite Element Method 

Perhaps the earliest documented large deflection, shell analysis program available was MARC, 
Ref. [1]. This program has grown over the years to include elements not exclusively for shell of re- 
volution analysis, and has become a general system for large deflection, inelastic analysis of struc- 
tures, Ref. [2]. The program contains both isotropic and kinematic hardening theory, and is able to 
represent the material properties with linear hardening, and Ramberg-Osgood representations. The 
solution scheme utilizes the tangent modulus method and Gaussian elimination. The basic revolved 
elements used in this program are those due to Khojasteh-Bakht, Ref. [3]. A second program using 
these elements is due to Yaghmai, Ref. [4]. However, this was a small program (limited to 80 ele- 
ments ) containing isotropic hardening theory. 

Another shell of revolution, finite element code (IBLAD) was developed by Zudans, Ref. [5]. The 
capability also included isotropic and kinematic hardening theories, with the solution process utilizing 
the tangent modulus method and Gaussian elimination. The experiences gained in formulating the above 
program were then utilized in creating a more general shell analysis program, EPACA, Ref. [6]. This 
program has the capability of elastic-plastic-creep analysis of thin or thick three-dimensional shells 
made up of curved or flat shell elements. Continuous or discrete elastic foundations are also included 
in the analysis. The program contains both isotropic and kinematic hardening theory, and is able to 
represent the material properties with linear hardening, and Ramberg-Osgood representations. The 
solution scheme utilizes the tangent modulus method and Gaussian elimination. Heat transfer analysis 
capability is in a developmental stage for this program, as is a completely new code, DYPLAS, Ref. [7], 
which will include dynamic effects. 

Other finite element programs, tailored more exclusively to shells of revolution include those of 
Ref. [8,9]. The DYNAPLAS program, Ref. [8, 10], is a small (limited to 50 elements and 20 Fourier 
harmonics) but sophisticated program which includes dynamic effects. It utilizes a modified Novozhilov, 
Ref. [11], shell theory, and contains both isotropic hardening and a mechanical sublayer model, Ref. 

[12]. The solution process utilizes the “right hand side” (initial strain or stress) technique, and a 
restart capability is available for longer problems. Rings and spring foundations can also be included 
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in the analysis. The AXSHEL program, Ref. [9], is actually a subprogram of the PLANS system, Ref. 
[13]. This system was designed and built specifically for nonlinear analysis (plasticity and large de- 
flection) of structures, and includes, among others, special crack tip elements. The program uses 
kinematic hardening theory, and is able to represent the material properties with linear hardening, 
and Ramberg-Osgood representations. The solution schemes utilizes the right hand side technique and 
partitioning. 

Other finite element, inelastic codes, which mainly do not use shell of revolution elements, but 
which nevertheless deserve mention, include ANSYS, Ref. [14], and NEPSAP, Ref. [15]. 

Finite Difference Method 

Several special purpose finite difference inelastic analysis programs have been available for 
sometime, Ref. [16-18]. The most general and recent of these, however, is PETROS2, Ref. [19, 20 ]. 
This program analyzes arbitrary thin shells using a two-dimensional finite difference space mesh, and 
expanding time also in central differences. First order difference formulas are used for all expan- 
sions. The mechanical sublayer model is used with integrations through the thickness being performed 
by Gaussian quadrature. Dynamic effects are considered. 

A finite difference program recently being developed for inelastic shell of revolution analysis is 
BOSOR5, Ref. [21]. The numerical method used should actually be classified as variational finite- 
difference, Ref. [22, 23], because the finite-difference approximations are applied to the integrands in 
the energy integrals, rather than the differential equations of equilibrium. The program contains only 
isotropic hardening although both flow theory or deformation theory may be utilized. Secondary creep 
effects are also included. The solution scheme utilizes the tangent modulus method with Newton-Raph- 
son iterations being performed in each cycle. 

Numerical Integration 

The only program using the numerical integration procedure for inelastic shell analysis to date 
was NONLEP, Ref. [24]. The program utilized Sanders shell theory, Ref. [25], and contained only iso- 
tropic hardening. There was no capability either for branched shell analysis or the solution of the 
coupled nonlinear torsion problem. The solution procedure for the nonlinear problem was Newton itera- 
tion. The newly available numerical integration program, STARS-2P, is the subject of this report, 
and its capabilities were briefly outlined in the introduction. 

The comparative advantages and disadvantages of the three numerical approaches above were dis- 
cussed in Ref. [26 1 and will not be repeated here. It should be noted, however, that the previous com- 
ments as to idealization grids should recieve even stronger emphasis in inelastic analysis, due to the 
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incremental nature of most solutions. 


A final note on theory accuracy is also in order. The computational capability available for the 
inelastic analysis of shells has experienced a phenomenal growth recently, as exemplified by the pro- 
grams discussed above. Indeed, the level of analysis capability that has been obtained has actually 
outstripped the ability to describe and verify complex material behavior such as cyclic, time and tem- 
perature dependent plasticity. The user of such programs must now be doubly on guard in using them 
in areas where the incorporated material behavior laws have not been verified by experiment. 
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CHAPTER 1 


GENERAL NONLINEAR SHELL EQUATIONS 

A full derivation of the pertinent equations may be found in Ref. [26] and [27], For the axisym- 
metric small strain-moderate rotation case they are (see Fig. 1 for the geometry and sign conven- 
tions): 


Strain-Displacement Relations 

The nonlinear axisymmetric strain-displacement relations, for shells of revolution, are: 
€ «s = e s« + 

i 2 

^ <p<p ~~ © <p<p 2 ^ 8 

€ ct = i(w 8 + w|) 

^ (p 8 ~ y Qv ^ 8 ^ ip 

where 


e 00 - e 88() ~ £k 8 


y0«.=y0<» o -2Sk Se 

and 

e 8 0 O = ^ {v cos<p - w sin<p} 
- 1 r i 

e «”o'7 1 l v . «-- w l 
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k 8 = - — cu 9 coS(p} 
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with w 8 and given by 


w 8 - — (w, „ +v) 


w 0 = - — (usimp) 
r o 


(4) 


Stress Resultant -Strain Relations 


The stress -resultants are defined as (see Fig. 2 and 3 for sign convention): 



+z 



8 






N,= ja,dt 


M, = 



N a9 - N Sp - 




( 5 ) 


The thermal resultants are defined by 


N Tfl = | 

N t ,= j 
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and the inelastic resultants are given by 

{ " J G " y «- 



E B (e p B + i , Bfl eg) {ll 

1-^8^ U) 


(6) 


(7) 


where the superscript, P, indicates plastic strains. 

Using the stress-strain relations fori an orthotropic body in a state of plane stress, and' using the 
centroidal surface as the reference surface of the shell wall, the stress- resultant strain expressions 
become: ! 

= K n [e 9 8 0 + v ev e VVQ ] — N T8 - N le , 

r , '\ 

N„ - K 22 [e 0B() + v a ae sfi(l J - N To - Ni„ 
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M b9 - M 9o — 2D 33 k B9 - M IaS 


Variations of the above equations may be obtained to account for different classes of smeared reinforce- 
ment, Ref. [26]. The definitions of the extensional stiffnesses K tl , K 22 , bending stiffnesses D u , D 22 , 
and in-plane shear stiffnesses K 33 , D 33 , are available, for a variety of cross sections in Ref. [26]. 

Equilibrium Equations 


The axisymmetric nonlinear equilibrium equations for shells of revolution are: 


1 

r o 
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(9) 
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r 2 
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( 10 ) 


The f* terms are the nonlinear contributions to the stiffness matrix (or they can be used, as shown, 
as effective geometric forces). The F, terms are the applied surface forces tangential and normal to 
the deformed shell surface, the f , are the forces in the undeformed coordinate system, and m„ and m 9 
are the distributed moments per unit surface area. It should be noted that the sixth equilibrium equa- 
tion cannot in general be satisfied with this shell theory because of the assumption that £ « r. How- 
ever, this equation is not employed in the subsequent solution of shell problems. , 

Boundary Conditions 


On a cp = constant edge, we must have one of each of the following boundary conditions specified. 


u or T w9 
v or N„ 

(U) 

w or J*=Q„ - r t fj 
or M„ 

with the following definitions for the “Kirchoff-type” effective stress resultants: 


T* 9 =N„ 9 - simp 
r o 

J® = Qe — N w9 W „ — N„ to 9 


(12) 
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Final Equations 

These above equations for the axisymmetric case may be manipulated and put in the form of eight 
simultaneous first order differential equations in y?, convenient for use in a numerical integration com- 
puter code, Ref. [26], For simple shell wall cross sections these are 


'«»,«> cos cp 
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(13n) 
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(14) 


For smeared reinforced wall cross sections more complicated sets of equations will arise, Ref. [26], 

In order to apply the incremental Lagrangian approach, the typical variable in Eq. (13, 14) may 
be defined as 

Yc w = Yc, + Ya (15) 

where Y c is the converged value of a stress or displacement in equilibrium with a specific value of 
load, and Y A is the increment necessary to accommodate an additional fraction of load. Since both 
Y c and Y Cj must satisfy Eq. (13, 14), we may substitute Eq. (15) into Eq. (13, 14) and subtract the 
form of Eq. (13, 14) containing only Y Ci to obtain the necessary differential equations which must be 
satisfied by Y A . Performing this operation and further simplifying by neglecting nonlinear terms in 
Y a , we can obtain the following incremental equations 
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In writing the above equation sets [Eq. (16, 17)] the subscript, A for incremental, has been omitted for 
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notational convenience. The R, are remainder terms which will be discussed in the section describing 
equilibrium corrections (Chapter 3). 

Equations (16) and (17) are the incremental equations which must be solved. Note that all the 
first order cp derivatives [Eq. (16a)-(16h)] are on natural boundary quantities, and that the equations 
contain previously accumulated stress and displacement quantities (subscript C ( ). 
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CHAPTER 2 


PLASTICITY ANALYSIS 

This section considers appropriate incremental plasticity relations to determine values of stress 
and plastic strain developed during the history of loading. In STARS-2P, two hardening rules, kinemat- 
ic hardening (see Fig. 5) and isotropic hardening (see Fig. 4), are available for initially isotropic, strain 
hardening materials. These are used in conjunction with the Von Mises yield condition. Drucker’s 
postulate is used to obtain the associated flow rules. For initially orthotropic materials, only kinematic 
hardening theory is available for strain hardening materials. It is used in conjunction with Hill’s ortho- 
tropic yield criterion. Ref. [28]. Ideally plastic behavior for initially isotropic or orthotropic materials 
may also be used. Additionally linear and nonlinear strain hardening capability is included in the pro- 
gram. Nonlinear hardening is achieved through use of the Ramberg-Osgood power law representation, 
Ref. [29], of the actual stress-strain data. A table of available plasticity options is presented in Fig. 

6, while Fig. 7 iUustrates how the input parameters determine which option is chosen. A detailed dis- 
cussion of how plastic effects are included, the matrix plasticity relations, kinematic and isotropic 
hardening, and how to best choose material constants foUows. 

Inclusion of Plastic Effects in STARS 

In STARS, plastic effects are incorporated using the “initial strain” concept, Ref. [30]. Plastic 
strains are interpreted to be initial strains present in the body in the same manner that thermal strains 
are treated. These initial strains are then applied as distributed “inelastic forces and moments.” This 
effective load technique differs from the “tangent modulus” approach in which the actual structural 
material properties are adjusted to account for plastic effects. For the initial strain approach the addi- 
tional assumption is made in the development of the governing differential equations [Eq. (13)-(14) of 
Chapter l] that these initial strains are independent of the displacements. The increments of plastic 
strains (initial strains) are then calculated using subsidiary plastic relations (to be presented later). 
These plastic strain increments are then used to calculate the inelastic forces and moments Nj and M t 
of Eq. (7). 

As discussed in Chapters 1 and 3, an incremental predictor procedure is used to calculate the 
displacement increments for the current step. It is assumed that the initial strain increments used in 
the current step are those calculated at the end of the previous step. At the end of the current step an 
“equilibrium check ” is made and any errors incurred in the prediction of plastic strains are adjusted 
through the application of the “equilibrium correction” in the differential equations, Eq. (16, 17). 

To determine the inelastic forces arid moments, the plastic strains are evaluated at each Runge- 
Kutta integration point at a specified number of layers through the thickness. For solid cross sections 
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a 1 Loading 



Fig. 4 Isotropic Hardening 



Fig. 5 Kinematic Hardening 
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Initial Material Behavior 
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Orthotropic 
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• Isotropic 
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Perfect Plasticity 


Note: Linear or nonlinear strain hardeninq is chosen on the 

. basis of input plastic parameters. r 



Fig. 6 Available Plasticity Options 
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Fig. 7 Available Hardening Options 
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a maximum of 20 layers (21 integration points) through the thickness may be specified and the number 
of layers must be even. Ten layers is sufficient for most problems. For sandwich cross sections the 
maximum number of layers that can be specified is eight (nine integration points) and both top and 
bottom face sheets (even if they are unequal in thickness) will contain this number of layers. The in- 
cremental plastic forces and moments are then obtained by numerically integrating the plastic strains 
through the thickness using Simpson’s rule. Ref. [3l]. 

Matrix Relations -Strain Hardening 

We can, for small strain increments, decompose the total strain increment {Ae T } into elastic 
{Ae 6 } and plastic {Ae p } components, as: 

{Ae T }={Ae 6 }+{A€ p } (18) 

The elastic strain increments are related to the stress increments {Acr} by: 

{Ae e }= [Ej-^Aa} (19) 

where [E]" 1 is an array whose elements are combinations of elastic material constants. A linear in- 
cremental constitutive relation between plastic strains and stresses can be written for all of the 
popularly used flow theories of plasticity. This relationship can be expressed as: 

{A€ p }= [C]{Aor} (20) 

Therefore, substituting Eq. (19) and (20) into (18) we can write: 

{Aor}=[R]- 1 {Ae T } (21) 

where [r]= [e]‘ 1 + [C]. 

The elements of the matrix [C] depend on the current state of stress, the yield condition, flow 
rule, and hardening law. In STARS the yield condition is based on Hill’s yield criterion for an initially 
orthotropic material, Ref. [28], which reduces to the familiar Von Mises yield criterion for initially 
isotropic materials. The associated flow rule is based on Drucker’s postulate, Ref. [32], Two harden 
ing laws are available. Prager- Ziegler kinematic hardening theory, Ref. [33, 34], is available for 
initially isotropic and initially orthotropic materials under cyclic loading. Isotropic hardening theory, 
Ref. [35], is available for only initially isotropic materials and is not suggested for cyclic loading 
situations because of its inability to predict the Bauschinger effect. 

Both linear and nonlinear strain hardening rules may be used with either hardening law. Input 
parameters determine which law is chosen (see Fig. 7). To minimize the input requirements for non- 
linear hardening, a Ramberg-Osgood, Ref. [29], representation of the stress-strain data is used: 
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a 3 a / g V' 1 
"E + 7E V 0 . 7 / 


(22) 


Thus for this representation of the stress -strain law, for an initially isotropic material two addi- 
tional material parameters, n and ct 0<7 , are required for the plastic analysis. For a Ramberg-Osgood 
representation of an initially orthotropic material these two quantities are required for each principal 
material direction and for each of the shear components. For a plane stress problem, as encountered 
in STARS, this means that three sets of stress -strain data must be specified (two for the normal com- 
ponents a v ,cr 9) and one for the shear component t„ s ). 


Kinematic Hardening 

Kinematic hardening theory as presented by Prager and Ziegler, Ref. [33, 34], predicts a rigid 
body translation of the yield surface (see Fig. 5) without any expansion or rotation. It predicts an ideal 
Bauschinger effect and was originally derived for use with linear strain hardening. In fact, one of the 
unresolved questions in the use of the kinematic hardening theory of plasticity is the definition of the 
hardening coefficient, c, under multiaxial stress states. This is especially true for problems involv- 
ing nonlinear hardening under cyclic loading. Several definitions have been proposed, Ref. [26, 27], 
with moderate success for initially isotropic behavior. An extension of this definition for initially 
orthotropic behavior based on Hill’s yield criterion has been recently proposed, Ref. [38]. 


Before discussing the hardening coefficient, let us first rewrite the basic equations required for 

a plastic analysis using kinematic hardening. 

\ 

First we need a yield condition, i. e. , , a rule to determine when plastic flow is initiated. In par- 
ticular, let us use Hill’s yield criterion for orthotropic materials: 

\ 

. f = F(a y -a z ) 2 + G(a I -o : x ) 2 +H(cr I -u 7 ) z + 2LT 2 „+2MT| x +2Nf 2 xy =l (23) 

\ 

where 

' \ ■ 

1 11 /'. 

3^ = g+ h ^ = h+f ^=f + g • 

2L = ^ 2M = J§- 2N-;j3- 

V 

Here - a if , with a 1( being the components of the shift in the yield surface, and X, Y, Z beingthe 

yield stresses in tension in the principal directions of orthotropy, and R, S, T being the yield stresses 
in shear with respect to the principal axes. For initially isotropic materials this reduces to the fa- 
miliar Von Mises yield condition: 

* =(^f ^ *■«- * 2 o (24) 
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with cr 0 equal to the initial yield stress in tension. 

Stability requirements for the use of Hill’s form of the yield surface, i. e. , insuring that it re- 
main a “closed ellipsoidal” surface, require that the following criterion must also be satisfied by the 
yield stresses, Ref. [39]; 

FuF 22 -Ff 2 = 0 (25) 


where 

^ 1 „ 1 _ 1/1 1 1 \ 

Fu - x2 Fs-jj-^ f j?-- 2 (x 5 + y 2 " Z V 

If Eq. (25) is not satisfied by the yield stresses then some other yield criterion must be chosen 
(see Ref. [39] for example) and implemented in the code. 

The associated flow rule used in STARS-2P is based on Drucker’s postulate: 

dcf, = d\ ' (26) 

8<r„ 

where dx is a positive scalor quantity. For kinematic hardening it is given by: 



(27) 


and c is the hardening coefficient to be defined. The incremental shift in the yield surface da,, is as- 
sumed to be along a radial line through the center of the yield surface as postulated by Ziegler, Ref. 
[34]: 

d a u = dp.(cTij - <*ij) (28) 


and 


dp = 


« j 

-r— dor i j 


, . 8f 


(29) 


In evaluating the terms of Eq. (26)-(29) for subsets of the nine dimensional stress space, all 
stress terms are retained in the yield function. Symmetry is applied, as well as setting stresses iden- 
tically equal to zero, after the expressions for dx and dp have been obtained. This is termed complete 
kinematic hardening. 


For plane stress kinematic hardening the appropriate terms for the relationship between plastic 
strain and stress increments given by Eq. (20) are 
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(30) 




■" Mf M,M 2 M t M 3 ‘ 


^=5 


M 2 Mj 


Mf 


M,M, 


_M 3 M! M 3 M 2 Ml J 


D = c^Ml + Ml+^+M^ 


where 


M x = 2(G + H)o il , - 2H6& 


M 3 = 4Nt v g 


M 2 = 2(F + H)a e -2H a v 
M 4 = — 2F(7|9 - 2Gd„ 


and the expression used to obtain the Aa u from Eq. (28, 29) is 


'AaA 


(o„ - 


M !Ao„ +M 2 Ao fl +M 3 At bS 

<0g — 0t e 

2 

■ Aa J 


v* - 


There still remains the problem of determining c, the hardening coefficient. If we evaluate c for 
a uniaxial state of stress and restrict ourselves to initially isotropic materials, using Eq. (26), (27), 
and (24), we get: 


1 = 3 de p 
c 2 do 


(31) 


If we generalize to a multiaxial state of stress, we may consider o and e p to be effective quan- 
tities o and € p , Ref. [36], where for an initially isotropic material 


+ +3t ~ +3t2 « + 3 t* (32) 

and 

d? p = V(2/3)de p de p 

Note the difference in the definition of the “effective stress” d and the yield function “l” of Eq. 

(24). 

If we use a Ramberg-Osgood representation of the stress -strain curve and consider the nonlinear 
term to represent the plastic strain contribution, then: 
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(33) 


1 3 de p _3^3n 

c 2 do 


n-1 


'3n 

a 

07E 

<70.7 


or for linear strain hardening 


1 _ 3 d? p 
c 2 da 


e t /e _ 


(34) 


where E T is the tangent modulus. 

Various methods of choosing an appropriate value for E T /E from a stress-strain curve have been 
suggested (see Ref. [40] for example). The use of linear strain hardening can lead to substantial un- 
derprediction of plastic strains in the range of strain where the material does not have a linear stress - 
plastic strain relationship. In turn this can lead to serious discrepancies if the response to cyclic 
loading is desired. 

To determine the values of n and a 0 . 7 that best fit the actual stress -strain data, the method sug- 
gested by Ramberg and Osgood can be used if the strain range is sufficiently small. For this case 
u lj log(17/7) 

Iog(ao.7/ CO. 85) 

The quantities a Q-7 and a Q-85 are the stresses at which the curve has secant moduli of 0.7E and 

0. 85E, respectively. If the strain of interest is sufficiently large so that the parameters as deter- 
mined by the preceding process do not fit the curve well, then a power law representation to fit the ac- 
tual data can be used 

€ = | + 0a n (35) 

Now having determined P and n to “best fit” the experimental data, the value of n input is that 
calculated, and the value of ao. 7 input is gotten by equating the nonlinear terms of Eq. (35) and (33), 

1. e., 


1 

n-1 



For cyclic loading the problem of how to determine c for load reversals subsequent to the initial 
loading is even more vexing expecially for nonlinear hardening. For linear hardening a procedure for 
selecting the value of c and new yield stresses for succeeding cycles is suggested in Ref. [40]. For 
nonlinear hardening it is desirable to reproduce the actual hysteretic stress -strain curves obtained 


25 



from test data. Morrow, Ref. [41], details extensive cyclic tests on isotropic metals to determine the 
cyclic and hysteresis stress -strain curves. He has found that both the cyclic and hysteresis curves 
may be well represented by a power law similar to the Ramberg -Osgood curves. Based on his find- 
ings, we have incorporated into STARS-2P a hardening coefficient for load reversals following the ini- 
tial monotonic loading defined as: 


c 2 da 


[3n 

a 

|_7E 

.2o‘ 0# 7 J 


where we now define a to be 


a 2 = 2 ~ - ?z * ) 2 + (— ~~) Z + 3t* 2 + 3r *f + 3 t 


*2 


(37) 


(38) 


where ofj =a l} - a li; and a u is the last value of stress at the end of the previous loading range. 


A generalization of the definition of c for initially orthotropic materials based on Hill’s yield 
criterion is presented in Ref. [38]. The basic highlights will be reproduced here. 

We define c, for plane stress, to be: 


1 1 do v ° v 1 3 1 dT„g fve 

c _ c„ 2/ + Cg 2/ + c v9 2/ 


where 


\ 


(G + H) 2 d<r„ 

C * " [(G + H) 2 + H 2 + G 2 J di^ 

(F + H) 2 dcTg 

C ® IIfThPTPTf 2 ] def 


c „ e =2 


d T v g 

dv^ff 


\ 

\ 



(39) 


This form for c reduces to the individual uniaxial stress-strain laws and the initially isotropic case 
given in Eq. (31). 


For a linear strain hardening approximation, we use: 

de% 1 ( 1 ~ e t „/ e ») del* 1 (1~ E T 9 / E g) 

dcr v ~ E v (E t#j /E„) ’ do 9 ~ E s (E Tg /E 9 ) 

d 1 (1- G T„/G g9 ) 

d G b9 


/ 


( 40 ) 


where E t#j , E Te and G T<;g are specified by the user for each component direction. / 
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For a Ramberg-Osgood form of the uniaxial stress-strain data we use: 


3n<, 


7E„ 


3na 


7E S 

°O.70 






1 <P& 


3n w9 


7 G vB 

r 0.7 vB 




where we have defined the effective stress a (for a full three-dimensional stress space) to be: 
3 


. 


2(F + G + H) 


[F(n y - crj 2 + G(o K - <r x ) 2 + H(n x - a T f + 2Lt^ + 2Mr l, + 2Nt 2 ] 


( 41 ) 


(42) 


for initially orthotropic materials. For plane stress, x corresponds to <p, y to 8, and all stresses 
with subscript z are zero. . • 

Now for plane stress, n„, a 0 . 7 v , fy, Oo. 7 S > n v e> and T o. 7 o9 are specified by the user for each com- 
ponent direction. 

For cyclic loading we specify <r to be: 

2 (f Tg + H ~ ) [F( < - °* )Z+G (^ -^ )2+ - <*>* +2Lt* s + 2Mt^ + 2Nt* 2 1 (43) 

where of s = <r u - 

Data for orthotropic cyclic loading is seriously lacking. If the representation given by Eqs. (39, 
41, 43) is deemed valid, then a 0>7ij should be set to twice the measured monotonic value for load re- 
versals after the initial loading range. Additionally at the end of each half cycle of loading new plastic 
material properties may be input so that each hysteresis loop may be accurately represented. 

At times, due to nonproportional loading or other reasons, local unloading occurs even though 
the applied load is increasing. The unloading criterion, which is checked in the program for every 
load increment, is given by 

da u so for loading dr neutral loading 

8o ij 


dn u <0 for unloading 


For plane stress kinematic hardening, this is given by: 
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(o 0 - |a e )dCT tJ + (a B - a e + 3f v6 dT 0 g < 0 


(44) 


for initially isotropic materials, and by 

[(G + H )o v - Ha e ]da 0 + [(F + H)ff s - Hffjdo, + 2Nr e , 9 dr o6 < 0 (45) 

for initially orthotropic materials. The values of dCT fJ used are obtained from the “elastic” stress- 
strain relations for that increment. They are the actual dcr^ if unloading is detected, and are neglected if 
continued loading occurs. For continued loading the actual d(T u are determined from the plasticity re- 
lations. If unloading is detected at a point, then all further stress and strain increments are elastic 
until reloading is detected using the appropriate yield criterion. 

Isotropic Hardening 

Isotropic hardening theory assumes that during plastic flow the loading surface expands uni- 
formly about the origin in stress space, maintaining the same shape, center location, and orientation 
as the yield surface (see Fig. 4). This representation of work hardening does not account for the 
Bauschinger effect exhibited by most structural materials. In fact, this theory provides that because 
of work hardening the material will exhibit an increase in the compressive yield stress equal to the in- 
crease in the tensile yield stress. For monotonic loading conditions, however, isotropic hardening. is 
satisfactory and is commonly used. In STARS-2P it may only be used with initially isotropic 
materials. 

The yield condition is the Von Mises yield condition given by Eq. (24) with the a,, equal to zero. 
Now ff 0 increases in magnitude to provide for expansion of the yield surface. For isotropic hardening, 
a 0 may be defined by Eq. (24) to be the "effective stress, a.” The flow rule is again given by 
Drucker’s postulate [see Eq. (26)]. 

For isotropic hardening we can show: 

dX = i^p (46) 

where H' = da/de e or the slope of the effective stress effective plastic strain curve. Also: 

. — (<r® — 2 xj$ )d(Tg + (oe — 20 g)d (T$ + ffdTv e 

da= = 

<x 


for plane stress. Here cr 2 = of -cr„o- 9 + of +3r 2 e . We can then write the elements of [c] as: 



M\ M,Mg M,M 3 
MjMg M| M 2 M 3 
M 1 M 3 M 2 M 3 Mf 


(47) 
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and 


D = <t 2 h' 

where 

yii~a 0 ~\o g , M 2 = ct s - \o v , M 3 = 3t 0S ’ 

These relations may be used with linear or nonlinear strain hardening laws as described in the 
section on kinematic hardening. Again the input consists of the value E x /E for linear strain hardening 
or n and cr 0-7 for a Ram berg- Os good representation of nonlinear strain hardening. It is not recom- 
mended that isotropic hardening be used for problems involving cyclic loading with stress reversals in 
the plastic range. If one wanted to do so, however, the plastic material properties may be changed at 
the end of each half cycle of loading. 

Local unloading is detected in a manner similar to that described in the section on kinematic . 
hardening, i. e. , 

8f , , 

- — da,,<0 for unloading (48) 

9o ij 

For plane stress isotropic hardening, this is given by 

(<r„ - i<r g ) do„ + (a„ - i<r v )da g + 3r vg dT v g <Q (49) 

where the values of da„ used are the potential elastic stress increments for that step. Again they are 
the actual da,, if unloading is detected and meaningless if continued loading occurs. This unloading 
criterion is identical to that used by Yamada et al. , Ref. [42]. 

Matrix Relations— Perfect Plasticity 

The treatment of multiaxial elastic ideally plastic behavior requires that the following conditions 
be satisfied: 

• The stress increment vector must be tangent to the loading surface. 

• The plastic strain increment vector must be normal to the loading surface, where the loading 

surface is the representation in stress space of the initial yield function or the subsequent 
yield function after some plastic deformation has occurred. 

If f(o„) represents the yield surface, the first condition can be expressed analytically as: 

~-da t , = 0 (50) 

oa„ 

and provides a linear relationship among the components of stress increment. Thus, one of the com- 
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ponents may be expressed in terms of the others. In matrix form, this can be written as: 

{Ao}=[E]]Acr} . (51) 

where {Ag } represents the independent stress components. 

The normality condition provides a linear relation among the various components of the plastic 
strain increment 

de y = dX (52) 

9cr,j 

This condition is derived from the flow rule and provides a linear relationship in which each of 
the components of plastic strain increment can be written in terms of any one component. This re- 
lationship may be represented in the following form: 

{Ae p }=[E]{Ae p } (53) 

where {Ae p } is the independent plastic strain increment. 

It is apparent from Eqs. (51) and (53) that the independent increments of stress and plastic strain 
can be combined and written in the components of a vector, {Aw} (see Ref. [36]) so that Eqs. (51, 53) 
can be written, respectively, as: 

{Ao}=[e]{Aw} (54) 

{*e P } = [E]{Aw} ' (55) 

Combining the above equations with Eqs. (18) and (19), we can form the following relation for the 
independent quantities : 

{Aw}= [E*].' 1 {Ae r } (56) 

where [E*] = [e]" 1 [e]+ [e]. 

These relations are used with Hill’s yield criterion [Eq. (18)] for initially orthotropic behavior 
which reduces to the Von Mises yield criterion for initially isotropic materials. Ideally plastic be- 
havior is specified in STARS- 2P by zero n and g 0i7 on the material property cards. For initially 
orthotropic materials either all components of the stress-strain law must be ideally plastic or none 
must be. 
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CHAPTER 3 


NUMERICAL SOLUTION PROCEDURE 


The matrix solution procedure utilized in the STARS-2P program is basically that of Ref. [26], 
and thus will only be briefly outlined here. More detail will be expended in discussing the specialized 
requirements created by the plasticity analysis; namely the plastic incremental load vector, and equi- 
librium corrections in an incremental analysis. 

The basic steps in the STARS numerical solution procedure can be outlined as follows: 


(i) For a numerical integration procedure, a segment length must be chosen. This is defined as 
the length along a shell meridian where numerical integration can proceed without encountering any 
small differences of large numbers in the subsequent matrix formulation. The differential equations, 
Eq. (16-17) are first integrated by the Runge-Kutta procedure, for each shell segment. Eight initial 
vectors are used on the homogeneous equations (Eq. (16-17)) with f t =N T1 =M T1 =N U =M n =Rj = 0), where- 
in each natural boundary quantity is set alternately to unity with all others being zero. A final integra- 
tion with all zero initial conditions is performed on Eq. (16-17), but now including all the load terms. 

(ii) The above integrations give rise to the following influence coefficient matrix equations be- 
tween the incremental shell forces, F, displacements, A, and loading, i: 


{F(j)} = [jFT]{f(j)} = [JFTK*! : 3C 2 : SC 3 ] 



{ 6(i) 

{A(j)} = [jDTK6(j)} = [JDT][y ! : : %] j f (i) 


(57) 


where the [JFT] and [JDT] are coordinate transformation matrices from the local (dashed) system to 
the global (solid) system (see Fig. 8), and the SC l; < y i matrices are obtained from the integration. Note 
that the effects of geometric nonlinearity are included in the 9C t , y t terms. The loading, £, consists of 
the applied distributed loads, thermal loads, and inelastic loads as defined by the inelastic stress re- 
sultants given in Eq. (7). The increments of plastic strain used to determine the quantities of Eq. (7) 
are taken to be those computed at the end of the current load step. These are used to predict the dis- 
placements at the end of the next load step. The equations above can be transformed into stiffness ma- 
trix form: 

{F(i)} = [lFT][-y 2 ]- l ([jDT] T {A(j)} - [<y 1 ][EDT] T {A(i)} - fo ,]{*}) 

{F(j)} = [ JFT][9C !][lDT]{A(i)} + [jFT][9C 2 ][-y 2 ]- 1 ([JDT] T {A(j)} 
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- to JmT] T {A(i)} - [%]{£» + [jftKsejM 


(58) 


(iii) Having thus defined the shell segment stiffness and load matrices from the integration pro- 
cedure, the program needs to know whether there is enough basic region capacity (1 segment = 1 region) 
to idealize the structure, or whether multi-segment regions are to be used. If the former is true, pro- 
ceed to step (iv). 

If multi-segment regions are to be employed, a matrix reduction procedure (Guyan, Ref. [43]) must 
be used. The segments are first grouped by order, and then the overall region stiffness and load ma- 
trices are reduced to matrices relating only the region end points. 

(iv) After completion of the region matrix manipulations, the direct stiffness matrix approach is 
used for the over-all problem solution. The results of such a solution procedure are the displacement 
vectors {a}, for all the region joints in the shell. Utilizing the matrices obtained during Guyan reduc- 
tion for each region, the displacement vector is expanded until it includes the displacements for all the 
segment ends throughout the shell structure. 

(v) Utilizing the displacement solution for each segment in Eq. (58) we can obtain the forces {f}, 
at every joint. Since the segment grid is allowed to be rather coarse, and is chosen on the basis of 
geometric parameters rather than load or expected stress distributions, a final numerical integration 
on the complete Eq. (16) and (17) is performed at this stage to obtain a full stress and displacement in- 
crement distribution within all the shell segments. The initial vector for this final integration is the 
{F(i)} and { A(i)} obtained from the over-all matrix solution. An automatic checking procedure is avail- 
able at this stage in the solution since the {F(i)}, { A(i)} obtained from the over-all matrix solution and 
the (F(j)}, (A(j)} obtained from the final integration must satisfy the appropriate joint equilibrium and 
continuity equations. 

Special Considerations Relating to Plasticity Analysis 

Since the current program uses an incremental formulation, both total and incremental values of 
the necessary displacement and stress resultant variables in Eq. (16-17) must be kept. Thus as each 
load increment is completed the accumulated values are updated using Eq. (15), and thence utilized as 
the Y Cl values in the next incremental pass through Eq. (16-17). Special calculations are also made due 
to the plasticity analysis as follows: 

(i) At the end of each incremental step in the analysis, inelastic loads as defined by the inelastic 
stress resultants given in Eqs. (7) are calculated if the cross section is in the plastic range. The in- 
crements of plastic strain used to determine the necessary quantities of Eq. (7) are taken to be those 
computed at the end of the current load step. 
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Three basic shell wall cross sections are incorporated into the STARS-2P program. These are a 
solid cross section, and a sandwich cross section with equal, and unequal, thickness face sheets. Plas- 
tic effects are included by evaluating the plastic strains at prescribed points through the thickness, and 
thence obtaining the “inelastic forces and moments” [Eq. (7)] by integrating the plastic strain through 
the thickness using Simpson’s rule. Up to 20 layers may be specified for the solid cross section. For 
sandwich sections, up to eight layers in each face sheet are allowed. The number of layers may be 
varied from segment to segment. These inelastic forces and moments are then incorporated into the 
next increment load vector, {f}, in Eq. (58), in the same manner as thermal forces and moments. 

For smeared rings and stringers, the stresses and strains are evaluated only at the centroid of 
the stiffeners in the direction normal to their cross section. The centroidal offset distance is specified 
by the user. For rings, inelastic forces and moments with respect to the shell reference surface are 
obtained as: 

N Ifl =EA H €e M IS = EA R z c e| > (59) 

These are summed with the plastic resultants obtained for the basic shell wall cross sections. To cal- 
culate the plastic strains in the rings, uniaxial stress-strain relations are employed. Smeared string- 
ers are treated in a similar manner except that the meridional strain at the stringer centroid is calcu- 
lated instead of the circumferential strain. 

For the plastic analysis of waffle and isogrid reinforcements, the two dimensional state of strain 
in the shell is transformed into the directions of the reinforcements, whose centroids are offset a pre- 
scribed amount from the shell middle surface (the Kirchoff-Love assumption that normals remain 
straight and normal is assumed to hold). Uniaxial stress-strain relations are then used to calculate 
plastic strains and stresses in the stiffeners. “Inelastic forces” are then calculated in each reinforce- 
ment direction and are back-transformed into the shell coordinates and added to those found for the ba- 
sic shell wall cross section. For all smeared reinforcement cases the inelastic forces and moments 
are incorporated into the segment load vectors, {f}, as noted previously. 

A detailed stress analysis of isogrid, waffle, or other smeared sections is not obtained using the 
previous approximations. However, general load paths and the redistribution of forces may be obtained 
from the analysis. In addition, critical regions requiring a more refined inelastic analysis will be in- 
dicated. If one is interested in the local behavior in the vicinity of stringers or waffle or isogrid con- 
struction, then a nonaxisymmetric analysis must be employed. 

In instances where detailed stress analysis is required, or the ring spacing is too large, the use 
of a “smeared” ring is inadequate. Discrete rings with nine possible cross sections are included in 
STARS-2P for such cases. These thin ring elements are based on curved beam theory, Ref. [44], com- 
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patible with the shell theory. The details of the linear elastic theory and its incorporation into STARS 
are presented in Ref. [26]. Transverse shear strains and cross sectional warping are neglected, how- 
ever, the shear center and centroid need not coincide. 

The plastic ring strains are calculated at discrete points on the ring cross section, and are in- 
tegrated over the area using Gauss -Legendre integration to determine inelastic forces and moments, 
e.g., 

N 1= f EfgdA Mr = J Eel’zdA (60) 

These are added to any thermal terms that are present, and incorporated as external loads acting on the 
ring. These loads are then used as part of the line load vector in the post Guyan reduction stiffness 
analysis. 

(ii) It has been noted, Refs. [8 and 45], that the use of an incremental formulation may require 
equilibrium corrections for larger nonlinearities. This type of correction is automatically provided in 
STARS-2P at every integration point. The accumulated stress and deformation state in the shell is ob- 
tained from Eq. (15). This state is then substituted into the total form of the differential equations 
[Eq. (13a -f) and (13h-j)], where all the quantities can be thought of as having the subscript C t . These 
equations are not identically satisfied because the inelastic strains used as effective forces to calcu- 
late the current displacements were assumed to be those computed at the end of the previous step. 

These predicted plastic strains are not exactly correct, leading to “remainders” to the identities [Eq. 
(13)]. These remainders are the R t terms which may be found in Eq. (16). They are “pseudo load” ad- 
ditives to the segment load vectors as a correction term to be utilized in the next step. In this fashion, 
total equilibrium and strain compatibility is maintained in the shell structure throughout the loading his- 
tory, and excessive drifting from the exact stress -strain curve is prevented. 
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CHAPTER 4 


NUMERICAL EXAMPLES 

To test the accuracy of the STARS -2P analysis several test problems were chosen and the results 
compared with available analytical and/or experimental data. 

Circular Plate Under Monotonic Loading 

First let us consider a circular plate under a monotonically increasing central load. This prob- 
lem was originally considered in Ref. [9] where experimental results up to a total load of 500 pounds 
were presented. The geometry of the plate is shown in Fig. 9, along with a load versus central deflec- 
tion curve. Also shown on the curve are the results from BOSOR5, Ref. [21], up to 1000 pounds. Ex- 
perimental load-deflection data is presented here up to a 1000 pound total load and the additional data 
was obtained from Ref. [46]. The material properties used in the analysis correspond to a Ramberg- 
Osgood representation of the stress-strain data, Ref. [9], given by: 

e = 10- 6 (0.0953ff+2.80xlO- 1 V’ 661 ) . 

The specimen was made of 2024-0 aluminum alloy. 

The results for STARS-2P compare extremely well with the experimental data up to the 1000 
pound value. The predicted deflections are always slightly greater than the experimental data, with 
the worst discrepancies of about 10% occurring in the vicinity of P = 250 pounds. This corresponds to a 
large portion of the plate being in the knee of the stress-strain curve. A load increment of 5. 5 pounds 
up to a maximum load of 1104 pounds was used. Another run using load steps of 11 pounds differed only 
slightly from the results presented, indicating a converged solution. The equilibrium correction term 
was used at every load step. Subsequent runs on a UNIVAC 1108 computer rather than an IBM 370, 
showed that the number of load steps could be halved again without significant accuracy changes on this 
machine. 

Buckling of a Spherical Cap Under Uniform External Pressure 

Next we consider a simply-supported isotropic spherical cap under a uniform external pressure. 
Wilkinson and Fulton considered the plastic buckling of such a cap in Ref. [47]. They assumed that the 
material obeyed a linear strain hardening law with a ratio of tangent modulus to Young’s modulus of 
0. 1 and yield stress to Young’s modulus of 0. 002. A shell geometry parameter X = 2[3(1 - v 2 )] 1/4 (H/h) I/z 
equal to 4.0 was chosen, where vis Poisson’s ratio, H, the maximum shell rise, and h, the shell thick- 
ness. Other geometric quantities are defined in Fig. 10. The load-deflection curve for this shell is 
also given in Fig. 10. Initial yield occurs at a load of 900 psi, and buckling occurs at a load of 2150 
psi which is 1.4% below the value of 2181 psi predicted by Wilkinson and Fulton. 
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Fig. 9 Load-Deflection for Centrally Loaded Circular Plate 
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Plastic Buckling of a Sandwich Spherical Cap 

Sharifi and Popov, Ref. [48], consider the buckling of a sandwich spherical cap under uniform ex- 
ternal pressure. They use a finite element analysis which accounts for the elastic core stiffness and 
transverse shear deformation, but the core is assumed to remain elastic and only the face sheets are 
allowed to go plastic. They compare their results with the theoretical linear (eigenvalue) analysis of 
Yao, Ref. [49], and some experimental work performed by Lin, Ref. [50], The latter comparison is 
inappropriate for the present purpose since it was observed that the cap buckled asymmetrically, and 
the results of Yao, Sharifi and Popov, and STARS -2P all assume an axisymmetric buckling mode. 

The geometry of the cap and its material properties are shown in Fig. 11. Also indicated in the 
figure are the buckling load calculated by Yao (P* = 1), the experimental buckling load (significantly lower 
than the axisymmetric mode), the load-deflection curve of Sharifi and Popov, and that calculated by 
STARS -2P. The deflection is that at a point located at r/a = 0.69, which is the point with the largest de- 
flection for this structure. The STARS analysis apex deformation results are also shown. Sharifi pre- 
dicts buckling at a load slightly below the linear (eigenvalue) analysis of Yao. The STARS -2P analysis 
predicts a buckling load of P* = 1.045 which is approximately 5.5% higher than that of Sharifi and Popov. 
This might be due to the fact that the core in STARS -2P is assumed to be a rigid spacer between the 
face sheets. Despite this difference in the model, the analytical results compare reasonably well. 

Plastic Buckling of an End-Loaded Clamped Circular Cylinder 

We next consider a circular cylindrical shell clamped at both ends and subjected to a uniform edge 
compression. Pertinent dimensions are given in Fig. 12, and the material is 6061-0 aluminum alloy. 
Material properties were derived from data obtained from the Materials Selector Handbook and verified 
in the Grumman Research Department Applied Mechanics Laboratory. The Ramberg -Osgood parame- 
ters used were E = 10 7 psi, u 0 .7 = 5050 psi, andn = 4.8. A proportional limit of 4300 psi was used. This 
cylinder was the subject of an experimental investigation performed at NASA/Langley Research Center. 
Figure 12 shows the experimental load-deflection curve up to initial buckling, Ref. [51], together with 
the one predicted by STARS-2P. The STARS program predicts a buckling load very close to that ob- 
tained experimentally, however, the predicted deflections are always larger than the experimental data. 

The load-deflection data for this structure is very dependent upon the shape of the stress-strain 
curve. Additionally it was found that stress-strain data for 6061-0 alloy can vary considerably depend- 
ing on material thickness and shape (i. e. , sheet or rod). Since the material data used was not obtained 
from the specimen tested, it is possible that the experimental and theoretical stress-strain curves dif- 
fered considerably. 
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Fig. 11 Elastic- Plastic Buckling of Sandwich Cap 
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Axial Load, Lb 



Fig. 12 Axial Load vs Cylinder End Deflection 
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Cyclic Behavior of a Centrally Loaded Circular Plate 


The final test problem considered is that of a circular plate centrally loaded by a circular punch, 
and resting on a circular outer edge support. The geometry of the plate and a diagram of the load his- 
tory is shown in Fig. 13. This plate was the subject of an experimental investigation performed at 
Grumman, Ref. [52]. The monotonic stress-strain curve for the 2024-0 aluminum alloy was presented 
in Ref. [9], The Ramberg -Osgood representation of the monotonic stress-strain curve used in 
STARS-2P for this problem is: 

e = 10' 6 (0. 0995a + 9. 634 x lO^V' 987 ) . 

Since no cyclic stress-strain tests for this material were available, it was assumed that the 
monotonic and cyclic stress-strain curves coincided. To get the size and shape of the actual hysteretic 
stress-strain loops the representation proposed by Morrow, Ref. [41], was employed (see Chapter 2 
for the calculation of the hardening coefficient). The comparison of the experimental and theoretical 
load versus plate central deflection curves is presented in Fig. 14 for the loading pattern shown in Fig. 
13. Kinematic hardening theory was used in the analysis. 

Considering that the cyclic stress -strain curve is unknown, the correlation is quite acceptable. 
There is an initial over -prediction of the deflection at the maximum load (P =1320 pounds). The elastic 
unloading curves are parallel and the shape and size of the hysteretic curves are quite similar if one 
allows for the initial overprediction. Near the minimum load (P = - 770 pounds) the theoretical strains 
increase more rapidly than do the experimental strains. This may demonstrate that the actual materi- 
al does not exhibit an ideal Bauschinger effect, and thus the predicted theoretical plastic strains are 
larger than the actual ones. Upon reloading to P = 1320 pounds, the theoretical results predict slightly 
larger deflections than those produced after the initial loading while the experiment shows slightly 
smaller deflections. The initial overprediction and consequent cycling between different strain ranges 
may account for this difference. Examination of the stresses and strains at the top and bottom sur- 
faces revealed that for this i-inch thick plate geometric nonlinearity was negligible. 
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Fig. 13 Circular Plate Under Cyclic Loading 
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APPENDIX 


CONVERSION OF U. S. CUSTOMARY UNITS TO SI UNITS 


The International System of Units (SI) was adopted by the Eleventh General Conference on Weights 
and Measures in 1960. Conversion factors for the units used in this report are given in the following 
table: 


Physical quantity 

U. S. Customary 
Unit 

Conversion 
factor (♦) 

SI Unit 
(**) 

Length 

in. 

0.0254 

meters 

Stress modulus 

ksi 

6. 895X10 6 

newtons/meter 2 

(N/m 2 ) 

Stress resultant 

lbf/in. 

175.1 

newtons/meter 

(N/m) 

Temperature change 

°F 

5/9 

Kelvin (K) 


♦Multiply value given in U. S. Customary Units by conversion factor to obtain equivalent value in SI Units. 
♦♦Prefixes to indicate multiple of units are as follows: 


Prefix 

Multiple 

giga (G) 

10 9 

mega (M) 

10 8 

kilo (k) 

10 3 

deci (d) 

10' 1 

centi (c) 

10’ 2 

milli (m) 

10' 3 


NASA-Langley, 1975 
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